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Abstract 

An important unresolved question in strong interaction physics concerns the pa- 
rameterization of power-suppressed long-distance effects to hard processes that 
do not admit an operator product expansion (OPE). Recently Bauer et al. have 
developed an effective field theory framework that allows one to formulate the 
problem of soft-collinear factorization in terms of fields and operators. We extend 
the formulation of soft-collinear effective theory, previously worked out to leading 
order, to second order in a power series in the inverse of the hard scale. We give 
the effective Lagrangian and the expansion of "currents" that produce collinear 
particles in heavy quark decay. This is the first step towards a theory of power 
corrections to hard processes where the OPE cannot be used. We apply this 
framework to heavy-to-light meson transition form factors at large recoil energy. 



1 Introduction 



An important unresolved question in strong interaction physics concerns the parameter- 
ization of power-suppressed long-distance effects to hard processes that do not admit an 
operator product expansion (OPE). For a large class of processes of this type the princi- 
pal difficulty arises from the presence of collinear modes, i.e. highly energetic, but nearly 
massless particles. 

In a series of papers |2j, f|, || Bauer et al. have developed an effective field the- 
ory framework that allows one for the first time to formulate the problem of soft-collinear 
factorization in terms of fields and operators, rather than the factorization properties of 
Feynman diagrams. Assuming that the field content of the effective theory correctly in- 
cludes all long-distance modes, the effective field theory formulation simplifies or renders 
more transparent the factorization proofs that form the conceptual basis of high-energy 
QCD phenomenology || [7]. The effective theory has been formulated only to leading order 
in an expansion in the hard scale up to now, as have been the diagrammatic factorization 
methods. 

In this paper we extend the formulation of soft-collinear effective theory to second order 
in a power series in the inverse of the hard scale. We give the effective Lagrangian and the 
expansion of "currents" that produce collinear particles in heavy quark decay. This is the 
first step towards a theory of power corrections to hard processes that do not admit an OPE. 
Previous approaches to this problem include the renormalon approach (reviewed in ||), 
where power corrections have been identified through the asymptotics of the perturbative 
series at leading power. This approach is severely limited since only those power corrections 
can be included that are related through renormalization to operators appearing already 
at leading power. The effective theory approach is free from this limitation. 

The outline of this paper is as follows: in Section || we set up the notation, power 
counting rules, fields and gauge transformations of the effective theory. In Sections |3| 
and f| we derive the effective Lagrangian of soft-collinear effective theory (SCET) and the 
effective heavy-to-light currents, respectively, to second order in the expansion parameter 
of the effective theory. We then discuss in Section [| power-suppressed effects that break the 
symmetries between heavy-to-light meson transition form factors at large recoil, extending 
previous work of two of us on perturbative symmetry-breaking corrections. We conclude 
in Section |. 

2 Power counting and fields 

We begin with detailing the power counting rules and fields from which the effective theory 
will be constructed. This discussion is a repetition of previous work [fj] , using our preferred 
notation, except for two aspects. Bauer et al. use a hybrid momentum-position space rep- 
resentation in their formulation of the effective theory, where collinear fields are labelled 
by the large part of the collinear momentum of the particles they create and destroy. The 
corresponding derivatives are replaced by "label operators". We find it more convenient, 
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Figure 1: Kinematics of heavy quark decay into a single cluster of collinear and ultrasoft 



especially for the extension to power corrections, to stay within the conventional position 
space representation, where the effective field theory can be formulated in more familiar 
field theory terminology. The second modification concerns the definition of gauge trans- 
formations and Wilson lines, and is necessary to extend the theory beyond leading order. 
The position space representation also requires that one expands ultrasoft fields in their 
position arguments. This will be explained in Section |3|. 

2.1 Kinematical definitions 

The effective theory is designed to expand a physical process in powers of a small ratio of 
scales, A. In perturbative computations the relevant scales are determined by specifying 
the external momenta of certain Green functions. In a real hadronic process these external 
momentum configurations are enforced by hadronic wavefunctions or the selection of spe- 
cific final states by the observer. The situation we are interested in here consists of a cluster 
of nearly collinear particles moving with large momentum of order m into a direction n_, 
such that the total invariant mass is of order (mA) 2 . This acquires a Lorentz-invariant 
meaning if there is a source for collinear particles. We may think for instance of a heavy 
quark of mass m at rest, which decays semileptonically into a set of partons with large 
energy and small invariant mass, as is relevant to exclusive decays, see Figure [I]. More 
complicated situations often arise in jet physics, where the final state can consist of several 
clusters moving into different well-separated directions. However, in this paper we restrict 
ourselves to the case of a single direction n_. 

We decompose the momentum p of a collinear particle as 



where are two light-like vectors, n 2 + = n 2 _ = with n + n_ = 2. The individual momen- 
tum components are assumed to have the following scaling behaviour, 



particles. 




(1) 



n + p ~ m, p± ~ mA, n^p ~ mA 2 , 



(2) 
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such that p 2 ~ m 2 A 2 as required. In the following we shall put m = 1 in scaling relations 
such as (||), understanding that the dimension of any quantity is restored by inserting the 
appropriate power of m. 

The invariant mass of the final state is not altered if we add to it a particle with 
momentum k whose components are all of order A 2 . The effective theory therefore also 
contains such ultras oft particles. The effective theory does not contain soft particles, whose 
momentum components are all of order A, as external states, since adding such a particle 
(with momentum k) to a cluster of collinear particles with momentum p implies an invariant 
mass (p + k) 2 ~ A in contradiction with our kinematical assumption. 

2.2 Fields 

We construct an effective field theory by introducing separate collinear and ultrasoft fields 
that destroy and create collinear and ultrasoft (anti-)particles, respectively. It follows 
that an ultrasoft field 4> us (x) varies significantly only over distances of order x ~ 1/A 2 , so 
derivatives acting on these fields scale as 

<9 M US ~ A 2 us . (3) 

Collinear fields have significant variations over n_x ~ 1, and x± ~ 1/A. Derivatives acting 
on collinear fields 4> c therefore scale as 

n + d<p c ~ 4> c , <9_L0 C ~ \(fi c , n_d<p c ~ A 2 c . (4) 

The effective theory is constructed to reproduce the Green functions of QCD under the 
kinematical assumption above as an expansion in A. With a large scale m present, quantum 
fluctuations with virtuality m 2 also contribute. These hard modes can be integrated out 
perturbatively. In practice, this amounts to Taylor-expanding the loop integrals in A, 
assuming that all loop momenta scale as m. The effective theory accounts for those loop 
momentum regions where this expansion breaks down, i.e. for those configurations where 
some propagators become on-shell as A — > 0. For the given kinematics this leads to the 
collinear, soft and ultrasoft modes defined above. Since soft modes do not appear as 
external states, they can be integrated out.[] This will be understood in the following, but 
since (with the exception of Section this paper is concerned with power corrections in 
A at tree level, the details of integrating out soft modes will not be important. 

In the case of QCD we introduce collinear and ultrasoft quark and gluon fields and assign 
a scaling rule to them. A collinear quark field has two large and two small components. 
The decomposition of the spinor is as follows: 

i/> c (x) = £(x) + r](x), £(x) = ^ c { x ), r](x) = t±^- ij) c {x), (5) 

J This is similar to non-relativistic effective theory where soft modes with momentum k ~ mv give rise 
to instantaneous effective interactions in the effective theory of potential heavy quarks and ultrasoft gluons 
and light quarks flC[] . 
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where (r/ T r/±)/4 are projection operators, and = rji+ri = 0. The scaling of these 
components is obtained from the projection of the QCD quark propagator. For £ fields 

(0\TZ(z)Z(y)\0) = ^ (0\TA(x)fov)\0)^ = J ^ ^_^ e -*.-y). (6) 

For collinear momentum p the right-hand side is of order A 2 . From this and the analogous 
equation for r\ with n_ and n + interchanged, we conclude 

£ ~ A, 77 ~ A 2 . (7) 

The small 77 field will later be integrated out. The corresponding argument for an ultrasoft 
quark field, which we denote by q, dictates 

q ~ A 3 (8) 

for all components of the spinor. Here k 2 ~ A 4 and d 4 k ~ A 8 for ultrasoft momentum has 
been used. 

The scaling of the collinear gluon field is determined by projecting the gluon propagator 
(in a covariant gauge with gauge fixing parameter a) 



(0|T^(xK(y)|0> = / 



d 4 p i 
(2^)4 p 2 + it 



-9- + (!-«)££ 
p A 



(9) 



on its various components. This requires the gluon fields to scale identically to a collinear 
momentum, 

n + A c ~ 1, A ±c ~ A, n_A c ~ A 2 . (10) 
For an ultrasoft gluon field, 

K s - A 2 (11) 

for all components. Comparing the scaling properties of derivatives with those of gluon 
fields we see that collinear and ultrasoft covariant derivatives have homogeneous scaling 
rules: the terms in a covariant collinear derivative iD£<f) c = (id 11 + gA^)<f) c acting on a 
collinear field all scale in the same manner. A corresponding statement holds for iD£ s (fi us = 
(id^ + gA^)^ acting on ultrasoft fields. 

In cases where collinear particles are produced in the decay of a heavy quark, we also 
need an effective heavy quark field. When collinear modes interact with a nearly on-shell 
heavy quark, the heavy quark will subsequently be off-shell by an amount of order 1. These 
off-shell modes are removed from the effective theory, and heavy quarks interact only with 
ultrasoft gluons. It is therefore appropriate to remove the large component mv of a heavy 
quark momentum p = mv + k just as in heavy quark effective theory by introducing the 
field 

h v (x) = l±ie imv - x Q(x) (12) 

which varies significantly only over distances x ~ 1/A 2 , since k is ultrasoft. The momentum 
space propagator of h v is l/(vk), so h v ~ A 3 just as the massless ultrasoft quark field. 
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We finally need to assign a scaling to the integration element d 4 x in each term of the 
effective action, which must reproduce the power counting for Feynman graphs in momen- 
tum space. Expressing each field through its Fourier transform, <j)i(x) = J d 4 pie~ lxPl <f)i(pi), 
we see that the integral over x in the action eliminates one momentum integral 

// n \ n „ /n— 1 \ n— 1 

d A x J] dA Pi exp(-ix 5» = / Ml rf V, (27r) 4 5( 4 )(p n + ]T ft ) (13) 
\i=l / i=l J \i=l / i=l 

in a term with n fields. The scaling of d 4 x is thus 1/A 4 if the eliminated integral / d 4 p n is 
over a collinear momentum, and 1/A 8 if it is ultrasoft. For products involving both collinear 
and ultrasoft fields, one must choose the eliminated momentum to be collinear and thus 
count d 4 x as 1/A 4 . This is because eliminating the ultrasoft integration in an integral such 
as / d 4 pi c d 4 p2c d 4 p3 US would no longer ensure that the momentum p 3 = —{p\ c + P2c) is 
ultrasoft. 



2.3 Gauge transformations 

The effective theory must have a remnant gauge symmetry from gauge functions U (x) that 
can themselves be classified as collinear or ultrasoft. 

It is clear that ultrasoft fields cannot transform under collinear gauge transformations, 
since this would turn the ultrasoft field into a collinear field. Under ultrasoft gauge trans- 
formations ultrasoft fields transform in the usual way. The collinear quark field £ is also 
multiplied by U under both collinear and ultrasoft gauge rotations. Since ultrasoft gluon 
fields can be added to collinear gluon fields without changing the character of the collinear 
field, ultrasoft gluons can be viewed as slowly varying background fields for the collinear 
ones. Then A c transforms covariantly under ultrasoft gauge transformations, while it trans- 
forms inhomogeneously under collinear gauge transformations, but with the derivative re- 
placed by the covariant derivative with respect to the background field. These properties 
are summarized by: 



collinear: 



ultrasoft: 



U c A c Ul + -U c 
9 

* A 



A, A r UL 



D us , U\ 



U us A us Ul + -U us [d,Ul 



q 



U us q, 



(14) 



where D^ s = <9 M — igA£ s is an ultrasoft covariant derivative. The heavy quark field h v has 
transformation properties identical to q. Note that the sum A c + A us transforms in the 
standard way under both types of gauge transformations. 

Under collinear transformations the components n + A c and A c ± acquire a component 
of order A 2 due to the ultrasoft gluon field in the transformation rule. Ultrasoft gauge 
transformations are homogeneous in A, i.e. all terms in the transformation law have the 
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same scaling in A, but operators such as iD± c £ transform inhomogeneously in A under 
ultrasoft transformations as well. A consequence of this is that operators that have a 
definite scaling with A are in general not covariant with respect to gauge transformations 
as defined above, but should be combined with higher-order terms in the A expansion into 
a gauge- covariant object. In the following, we will always keep gauge invariance manifest 
as far as possible, before we switch to the strict A expansion of the effective operators. 

2.4 Wilson lines 

The effective theory makes extensive use of various Wilson lines. We collect the relevant 
definitions and properties here. The basic Wilson line is defined as 

W(x) = Pexp (ig J dsn + A(x + sn + )^j , (15) 

where P denotes path-ordering defined as P[A(x)A(x + sn + )} = A(x)A(x + sn+), if s < 0. 
Here A = A c + A us denotes the sum of the collinear and ultrasoft gluon fields. The Wilson 
line is a unitary operator, WW^ = W^W = 1, and has the following useful properties 

in + DW = Win+d, (in+D + it)' 1 = W (in+d + ie)~ l Wl (16) 

Under general gauge transformations U(x), W(x) transforms as 

W(x) -> U(x)W(x)U\x -oon+). (17) 

The particular form of the transformation matrix at ri-X — > oo turns out to be irrelevant, 
since Wilson lines always appear in combinations where the transformation at infinite 
distance drops out. An example is the finite-distance Wilson line W(x)W*(x + sn + ), 
which corresponds to a line integral between x + sn + and x, and which transforms as 

W{x)W\x + sn + ) -> U{x)W{x)W\x + sn + )U\x + sn+). (18) 

It is therefore consistent to choose U(x — oo n+) = 1 and this choice will be adopted in the 
following. Under collinear and ultrasoft gauge transformations we then have 

W -»• U c W, W -> U w W. (19) 

We also introduce the ultrasoft Wilson line Z that is defined as in (fL5|) but with the full 
gauge field replaced by the ultrasoft one, A us . This Wilson line is invariant under collinear 
gauge transformations and transforms as 

Z^U US Z (20) 

under ultrasoft ones. A particularly important quantity is the product of Wilson lines, 
WZ' . It fulfills the equalities 

in + D WZ ] = WZ ] in + D w , (in+D + ie)- x Wtf = WZ ] (in + D us + xe)' 1 (21) 
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and transforms as 

WZ ] -> U c WZ\ WZ^ -> [/ us WZ f (22) 

under collinear and ultrasoft gauge transformations. Furthermore we specify the two lim- 
iting cases 

WZ^ -> 1 for A c -> 0, 

WZ f -> W c = Pexp fig y° dsn + 4(i + sn+)j for A us -> 0. (23) 

Note that the purely collinear Wilson line W c defined above does not have a simple be- 
haviour under gauge transformations as soon as ultrasoft gauge fields are taken into ac- 
count. 

The integration measure ds in the definition of the Wilson line counts as 1/ A 2 if the in- 
tegrand contains ultrasoft fields only, but it counts as 1, if the integrand contains a collinear 
field. This can be seen in the momentum space representation, where J ds corresponds to 
l/(n + p). It follows that all three Wilson lines, W, Z, W c count as order 1. In contrast to Z 
and W c , however, W is not homogeneous in A due to the occurrence of both, collinear and 
ultrasoft gauge fields. Due to the case-dependent counting rule for ds, W does not even 
have a useful expansion in A, where ultrasoft fields should appear as corrections. What 
we actually need below is the A expansion of the product WZ\ which is well-behaved. To 
construct the expansion we rewrite (|2l| ) as 

in + D c WZ ] = \WZ\ gn + A us ] , (24) 

where the derivative acts only on WZ', and solve this equation iteratively with the "bound- 
ary condition" WZ* — > 1 when A c — > 0. We obtain 

WZ\x) = W c (x) + W c (in + d + it)' 1 W\ [W c , gn + A us ] + C(A 4 ) 

= W c (x) (l-i J° ds Wl{x + sn + ) [W c , gn+A us ] (x + sn + )^j + C(A 4 ), (25) 

where we used that the inverse of in + d, here defined with a -Me-prescription, acts as 

1 ro 

x) = — i / ds 4>(x + sn+). (26) 



in + d + it 

Note that the commutator term in (p5[), [Wc, gn+Aus] = [W c — 1, gn+Aus], has at least one 
collinear field when the exponential is expanded, and therefore counts as a collinear field. 
The integration element ds then counts as 0(1), which is crucial to ensure that the second 
term in (|25| ) is A 2 suppressed. In a similar manner one sees that (1 — WZ>) is also a collinear 
field. On the other hand, the Wilson lines W c and WZ* without the unity subtracted do 
not qualify as collinear fields and expressions such as 

'"' ds[WZ^f us ](x + sn + ) (27) 



where f m (x) denotes a product of only ultrasoft fields, do not have a well-defined A scaling 
behaviour and must be avoided. 
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3 Effective Lagrangian 



In this section we construct the effective Lagrangian for collinear and ultrasoft quarks 
and gluons including all interactions suppressed by one or two powers of A. We shall 
first derive the Lagrangian in a form that is manifestly gauge-invariant, but where the 
operators do not have a homogeneous A scaling. Each operator generates a series in A 
and is counted by its largest term in the A expansion. In a second step each operator is 
expanded in a series of terms each of which has a definite A scaling. We then show that the 
effective Lagrangian obtained by tree level matching is not renormalized to any order in the 
coupling a s , provided a Poincare- invariant regularization procedure such as dimensional 
regularization can be chosen. This remarkable non-renormalization property extends to all 
orders in A. 

3.1 Collinear Lagrangian 

We begin with deriving the Lagrangian for collinear quarks interacting with collinear and 
ultrasoft gluons, neglecting the ultrasoft quark field for now. This will be corrected for in 
the next subsection. The QCD Lagrangian for massless quarks is 



where ijj c is a four- component spinor, assumed to describe a nearly on-shell particle with 
large momentum in the n_ direction. The covariant derivative is defined as D = d — igA, 
where A is the sum of the collinear and ultrasoft gluon fields. Inserting the decomposition 
(|) into the QCD Lagrangian, we obtain 



Here we indicated explicitly that the ie prescription of the original Lagrangian is attached 
to the transverse derivatives. Having this in mind, we will suppress it in the following. 

We proceed by integrating out the small field r\ from the path integral. Since the 
Lagrangian is quadratic this amounts to solving the classical equation of motion, which 
gives 



where we have used (0) and (|26|) to obtain the second equality. This is consistent with 
7] ~ A 2 as stated in (|7|). In writing the solution for rj we have defined the inverse of in + D 
with a +ie prescription. This prescription is not given by the QCD Lagrangian, and we 
could have chosen the opposite sign or a principal value prescription. The ambiguity in 
this prescription cancels out in physical quantities, but the prescription is needed to define 
individual terms that appear in the calculation. This is similar to the need to define a 
treatment of spurious singularities of gluon propagators in light-cone gauge. The functional 






(29) 




S 



determinant of the operator in + D, which appears upon integrating out 77, amounts to an 
irrelevant gauge-field independent constant and can be dropped. To see this, note that the 
determinant is gauge-invariant so that we can evaluate it in light-cone gauge n + A = 0, 
where in + D = in+d.f\ The result for the collinear quark Lagrangian is then simply given 
by inserting ( p0|) into (|2£f), and reads 



2 m + D 
£(x) in^D — £(x) + i 



ie 2 



ds 



X 



[x + sn A 



(31) 



where D = d + igA is the covariant derivative acting to the left. Note that the effective 
Lagrangian is formally not hermitian because of the +ie prescription on (in + D)~ 1 , which 
also implies that products of fields 4>{xi) ■ ■ ■ 4>(x n ) are "ordered" along the n + direction, 
n^x\ > . . . > n_x n . The hermitian conjugate of the Lagrangian corresponds to the — ie 
prescription. The non-hermiticity is of no consequence, since it is related to modes with 
n + p = 0, while collinear modes have n + p ~ 1. Concretely, the non-local quark-gluon 
vertices that follow from expanding the Wilson line in the gauge field individually contain 
factors l/(n + p + ie), but these unphysical poles cancel in the sum of all diagrams. 

The Lagrangian (|3lD is equivalent to the Lagrangian derived in || in the hybrid 
momentum-position space representation, when W is replaced by the collinear Wilson 
line W c and when the ultrasoft gluon field is neglected in ilp±. Using the power counting 
rules developed in Section || it is easy to see that this corresponds to the leading term in 
the expansion in A. The result ([31]) is the exact collinear Lagrangian to all orders in A. It is 
manifestly collinear and ultrasoft gauge-invariant, which follows from the transformation 
properties fllTf ) of W, or directly from the first line of (31). The expansion of (|3l| ) in A will 
be discussed below, after ultrasoft quarks have been included. 

The derivation of the Lagrangian did in fact not require the assumption that ip c or £ 
are collinear fields. It can also be viewed as integrating out two of the components of the 
quark spinor field in full QCD, with n_ and n + chosen arbitrarily. This makes it clear that 
the Lagrangian ( |3lD is equivalent to the quark Lagrangian of full QCD, including all hard 
and soft modes (and ultrasoft quarks), provided we consider £ as a field that destroys and 
creates all these modes, not only collinear ones. This peculiar result follows because the 
notion of a collinear particle has no Lorentz-invariant meaning in the absence of sources 
that create such particles. The collinear Lagrangian must therefore exactly reproduce the 
Green functions of full QCD (projected appropriately as in ([|)), when they are evaluated in 
a frame where all particles have large momentum in the n_ direction. This interpretation 
is also confirmed by the observation that the equation of motion that follows from ( |31~D 
and the expression for 77 is identical to what one obtains in a treatment of full QCD in the 
infinite momentum frame |Tl| or in light-cone quantization |12|, |13~ . In the following we 



2 In a general gauge the functional determinant reproduces the quark loop diagrams with only 77 fields 
and any number of external n + A gluon fields. The loop integral vanishes, because the propagator for the 



quark lines, (in + d 



implies that all poles lie on the same side of the real axis. 
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shall not pursue this interpretation further since the Lagrangian will eventually be coupled 
to sources of collinear particles and we are interested in performing an expansion in A, 
which implies that fields must be classified as collinear and ultrasoft. 



3.2 Including ultrasoft quarks 

We now incorporate the ultrasoft light quark field, q, into the effective Lagrangian. The 
QCD Lagrangian gives rise to the following tree-level interactions among collinear and 
ultrasoft fields: 

+ £gAq + vgAq + qgA£ + qgAv + qnp w q, (32) 

where iD ns = id+gA us . Note that terms such as £ ilp us q vanish by momentum conservation 
enforced by the d A x integral in the action. Collinear gauge invariance is not manifest in 
fl52p. We now integrate out the rj field as before and order the terms in an expansion in A. 
In this process we use the field equation for £, which is equivalent to a field redefinition such 
that the final Lagrangian will be manifestly gauge-invariant under the transformations (|T4]). 

Since the Lagrangian ( |32"D is still quadratic in the rj field, it can be integrated out 
straightforwardly. The rj field is now given by 

In the following (in+D)" 1 will always be understood with a +ie prescription. Inserting 
(0) into © results in 

1 rl 

c = c c + c us + £g4cq + qgfict, - qqA- — 7^^-gAi 

- 101. -^—^gA q - qgA n?±L (34) 

with C c given by (|3TD, and £ us = qilpus q the purely ultrasoft Lagrangian. Note that both 
C c and £ us are order one terms in the action, since the integration measure d 4 x counts as 
A -4 in / d 4 x C c but counts as A~ 8 in / d A x £ us , where the integrand involves only ultrasoft 
fields. The mixed terms in the Lagrangian are subleading contributions in A to the action, 
beginning at order A. 

The new terms at order A involving the ultrasoft quark field are 

= q + q[g^ c -gn + A c -^—iIpAi. (35) 

This can be simplified with the identity 

-J— gn+ A c( j> w = (i-Wtf)<f> w -^—(l-Wtf)in + D w <l> w , (36) 
m+D m+D 
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which follows from ([21]) and is easily verified by acting with in+D from the left. The key 
point to note about this identity is that the first term on the right-hand side is of order 1, 
but the second is of order A 2 due to the ultrasoft derivative. For this to be correct it is 
essential that (1 - Wtf) is a collinear field, which ensures that (in + D) 1 counts as 1 and 
not as 1/A 2 . Thus, neglecting terms of order A 3 , we can write 

(gA ±c - ilpx. gn+JU q = (gA ±c - ± (1 - WZ*)) q = (tD±WZ^ - iD ±us ) q. 

(37) 

The last term gives £Hp± us Q, which vanishes in the action by momentum conservation, so 
£W takes the final form 

£« = lHp L WZ ] q + qZW^iJP\X + 0(A 3 £ C ), (38) 

which is manifestly gauge- invariant. 

With a similar reasoning we can bring the terms of order A 2 in (j34|) into a manifestly 
gauge-invariant form. Consider first the term involving two ultrasoft quark fields, which 
using fl3~6f) can be approximated as 

1 ^ + " A .,qn -qgn + A c {l-WZ^)^q. (39) 



in+D 2 °' " " T v ' 2 

Writing gn + A c = in + D — in + D us , using the property (|2l|) of the Wilson line and momen- 
tum conservation, we see that this term is in fact of order A 4 and can be neglected. This 
leaves us with the remaining terms at order A 2 , 

= e f gn„ A c + ilpx. ^— gA ±c ) i± q + qi± ( gn „A c + g A ±c ^— iD ± ) £. (40) 
\ m + D ) 2 2 \ m + D J 

Consider the term with q to the left. Using momentum conservation gri-A c can be replaced 
by in^D. Then we can rewrite this term as 

q ZW^ in_D + g (1 - ZW^) in^D ^ + q gA ±c — ^- H?J±£. (41) 
2 2 m+D 2 

The point of introducing the Wilson lines is that we can now use the leading-order equation 
of motion for the £ field, 

in the second term. This would not be correct in a term with ZW' or 1 in place of 1 — ZW^. 
The field equation cannot be used when multiplied from the left by a field which is not 
collinear, because it has been obtained by functional differentiation of the action with 
respect to the collinear field £. Multiplying the field equation with an ultrasoft field would 
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project out the ultrasoft component of the operators in ( [42]) which is not consistent with 
their original appearance as collinear operators in the action. Now (|4l|) reads 

qZWUn^D^-qUl-ZW^i^-gfiJ r^^i^. (43) 

Z L J VYl-i-i-J L 

The square bracket equals 

'zwHlp ± - 1$ LUS ZWA + iip ±us (1 - ZW r ). (44) 



We can integrate the ultrasoft derivative by parts in the last term to see that it is higher- 
order in A. Here it is again crucial that (1 - ZW*) is a collinear field. On the contrary, 
the first two terms individually are not collinear fields and only the difference of the two 
ensures that (in+D)^ 1 in (f43|) does not count as 1/A 2 . This gives the final result 

qZWUn-D^t + q (ZWUl? ± -iJp^ZW*) -J-^ilpJ±i, (45) 

for one of the two terms of (0), which is now manifestly gauge- invariant. The other term 
can be manipulated in a similar way, where we must now use the leading-order equation 
of motion for £, obtained by taking the functional derivative of C c with respect to £,[] 

£m_S^± = -£itf> ± (in + t> - ie)^ i*tf> ± —, (46) 



2 s~*r± V -+~ -r± 2 > 

where we have introduced (in + D — ie) _1 = (in + d — ie)^ 1 W and 

<f)(x) (in+d — ie^j ds (f)(x + sn + ) . (47) 

Collecting all terms in the expansion we find our final result for the soft-collinear effec- 
tive Lagrangian, 

£ = i (in-D + iiPt-^—iipA + q z^ us q + lilpiWtfq + qZWH^ 
\ in + JJ J 2 

+ li± in _ D W Z^ q + £ % z#j_ — ^- (ilp ± Wtf - WZ*iIp ±w ) q (48) 
2 2 m + D v ' 

+ qZW ] in_D ^ + q (ZWU^ - iJp ±w ZW^) — HpJ^t + 0(X 3 C c ). 
2 \ / in + D 2 

The terms containing only £ and only q are leading in A, where one must remember that 
the integration element d A x in the effective action counts as A -8 for a purely ultrasoft term 
and as A -4 for all others. The remaining two terms in the first line are order A interactions 
and the remaining two lines are order A 2 . The effective Lagrangian is manifestly invariant 
under collinear and ultrasoft gauge transformations. The quark-gluon Lagrangian flUf ) is 
supplemented by the pure Yang-Mills Lagrangian, which at this first step of the A expansion 
retains its original form. 

3 For this we use / d 4 xg(x) (in+d + ie)~ 1 f(x) = — J d 4 xg(x) [in+d — ie)^ 1 f(x), which follows from ( |26| ) 
and @. 
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3.3 Expansion in A 



In deriving (f48|) we counted the operators by their smallest power of A. Each term in 
the effective Lagrangian (f|8|) gives rise to a tower of operators which have a homogeneous 
scaling behaviour, meaning that the insertion of such operators gives a contribution of 
order A n to the scattering matrix element, where n is a certain fixed number that can be 
determined from the power counting rules. In a typical application of the effective theory 
we calculate the scattering matrix element to some order in A, so we must determine the 
Feynman rules that follow from the Lagrangian order by order in A. The expansion in terms 
with a definite (homogeneous) scaling in A is also necessary to define effective theories with 
multiple scales in dimensional regularization 



TO 



We now discuss this further expansion of the effective Lagrangian (fig). We must first 
split iD± = iD± c + gA± us , the first term being of order A and the second of order A 2 , 
and expand the Wilson lines in A. The expansion of WZ* has already been discussed in 
Section |2j4]. The combination W(x)W^(x + sn + ) of Wilson lines along the finite distance 
from x + sn + to x, which appears in the Lagrangian, also has to be expanded in the ultrasoft 
gluon field. Contrary to W(x) itself, this combination can be expanded in A, since it stands 
between collinear fields so that s counts as order 1. We obtain with \s\ ~ 1 

W(x)W\x + sn + ) = W c (x)W^(x + sn + ) 

+ i [° dt W c {x) \WI gn+A us W c ] {x + tn + )W^(x + sn + ) + C(A 4 ). (49) 



We need this to expand the inverse differential operator as 



1 



1 



1 



in+D in + D c in+D 



-gn + A v 



1 



' in + D c 



£>(A 4 



(50) 



where iD c = id + gA c . This expansion holds if collinear fields appear both to the left and 
to the right of the operator. 

In addition we must expand the position arguments of ultrasoft fields in the transverse 
direction and the direction of n + , since the ultrasoft fields vary more slowly than collinear 
fields in these directions. Defining x„ = |(n + x)n„, the relevant expansion is 



<f>us{x) = 0us(z-) + x±d<fr us {x-) 



H — n-X 
2 



n + d<p ns + - x^Xu^d"^ 



[X- 



C(A 3 



(51) 



The collinear field multiplying this expansion varies over distances x± ~ 1/A and n_x ~ 1, 
hence the term with <9j_0 us is of relative order A and the terms with n + d<p ns and d±d±<j) us 
are of relative order A 2 , since all derivatives on ultrasoft fields scale as A 2 . In terms of 
momentum space Feynman rules this Taylor-expansion corresponds to the fact that in a 
collinear-ultrasoft vertex the momentum transfer by the ultrasoft lines to collinear lines 
is small in the transverse direction and the direction of n + compared to the momentum 
carried by collinear lines. To obtain a homogeneous expansion in A the momentum space 
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diagrams must be expanded in these small momentum components. As a consequence 
of this momentum is not conserved at the interaction vertices. In coordinate space this 
corresponds to a breaking of manifest translation invariance due to the Taylor-expansion 
of ultrasoft fields around an arbitrary point in the transverse and n + direction. Of course 
translation invariance is recovered order by order in A. 

Performing this "multipole" expansion up to order A 2 , the Lagrangian reads 

C = Cf + Cf + cf + cf + cf + £ us , (52) 

where the various terms are given by 



Cf = £\in_D + iip. 



1 




in + D c 

cf = £ (W± c — l Tr 94±us + Mius^tt- i#L C + [(x±d) (jn-An)]] % £, 
\ in+Dc in+JJc J l 

+ 1 [iIPlc—^- [(x ± d) gA ±us ] + [(x ± d) gA_ " : ' 




y in + D c in + D 



+ (n-x) i(n + d)(gn_A ns )] + \x^x\ [d,d u (gn-AJ)]^ £ (53) 
for the part of the Lagrangian without ultrasoft quark fields, and 



Cf = £il?XcW c q + qW}il?xcZ, 



Cf = £g4±mW c q + ^[in-D + il? ±c ^—il? ±c }W c q 

z \ in+Dc ) 

+ iiip± c W c [(x ± d)q] + u h.c:' (54) 

for the interactions between collinear and ultrasoft quarks. In the last line "+ h.c." refers to 
the terms of the form q[. . .]£. In these expressions, and in the remainder of this subsection, 
all ultrasoft fields are understood to be evaluated at position x_. For instance, in_D 
at point x equals m_<9 + gn_A c (x) + gn_A us (x^). The square bracket means that the 
derivative operates only inside the bracket, and that all ultrasoft fields or their derivatives 
are evaluated at x = X- after the derivatives are taken. On the other hand, in an expression 
such as Hp± c W c q, the field q is evaluated at x_ so that the transverse derivative effectively 
acts only on W c . Note that we cannot use the equation of motion for £ to drop the 

(2) 

second term in C^ q , since W c q to the right of the operator is not a collinear field. The 
terms in the Lagrangian without the ultrasoft quark field and without factors of x from 
the multipole expansion appear to be equivalent to the corresponding terms in the hybrid 



momentum-position space representation of the Lagrangian given in [13] on the basis of 



reparameterization invariance. The derivation presented here shows that (leaving aside 



14 



the terms generated by the multipole expansion and those involving the ultrasoft quark 
field which have not been considered in [14|) these are the only terms that appear in the 



effective Lagrangian at order A and A 2 . 

The individual terms 0- % > of the effective Lagrangian ([52]) are no longer invariant under 
ultrasoft and collinear gauge transformations. The transformations flI3|) mix terms of 
different order in the A expansion, so that only the sum Y^l=o £^ is gauge-invariant up to 
higher-order corrections. This is related to the inhomogeneous transformation of collinear 
fields under multipole-expanded ultrasoft gauge transformations, 



Ox) -> U ws (x)£(x) = U Wi {x-)t(x) + [(x x d)U us }(x„)Z(x) + ... 
A c (x) -> U us (x^A c (x)UUx_) + l(x ± d)U us }(x^A c (x)UUx-) 
+ U w (x-) A c (x) [(x ± d) uU(x-) + ■■■, 



(55) 



and similarly for W c (x). Furthermore, under collinear gauge transformations flU]) the 
collinear gluon projections (n + A c ) and A± c receive A suppressed contributions from ultra- 
soft gluon fields. The only gauge invariance that remains manifest for every term in the 
A expanded effective Lagrangian (|52|) is for ultrasoft gauge transformations £/ us (x_) that 
depend only on a;_. In this case collinear fields transform homogeneously with U us (x^). 
The same is true for partial derivatives [d±q], [(n+d) q] and ultrasoft gluon field compo- 
nents A± us , n + A us . Although the lack of manifest gauge invariance causes no problem for 
the calculation of on-shell Green functions in the effective theory, it would be desirable to 
cast the Lagrangian into a form which is manifestly gauge-invariant under the full class of 
gauge transformations (|I4|). This can be achieved by making use of the equations of motion 
for collinear quarks and gluons, which is equivalent to a redefinition of the corresponding 
fields in the effective Lagrangian. In the following we demonstrate how this can be done 
for the simpler case of abelian gauge fields. 

In the abelian theory ultrasoft and collinear gauge fields do not interact, and the 
collinear gluon field does not transform at all under ultrasoft gauge transformations. Us- 
ing the equation of motion for the collinear quark field is then sufficient to rewrite the 
effective Lagrangian in the desired form. The terms and do not contain ultrasoft 
gluon fields and already take their final form. The remaining order A corrections to the 
soft-collinear Lagrangian can be rewritten with the help of the commutator identity 



±c 



in + D, 



iip. 



±c 



This results in 



4" = f t 



(gx±A ns ), ijpx 



in + D c 



in + D c 



-ilpi 



in + D c 



+ [(x ± d) {gn„A w )\\^-t 



(56) 



(57) 



We can now use the leading-order equations of motion for £ and £ to obtain the final result 
for C 



(i) 



Cf = C (« gF™) 4± £ (abelian gauge field) , 



(5f 
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where F denotes the field-strength tensor. This interaction term is the light-front equivalent 
of the electric dipole interaction x ■ E of ultrasoft photons with non-relativistic electrons. 

Similar manipulations bring the second-order terms £f\ Cfjj into a form in which 
A us appears only in in^D or the field-strength tensor. To derive the new form of the 
Lagrangian, one has to keep track of the changes of the Lagrangian at order A 2 induced 
after every application of the equation of motion at order A. The details of this derivation 
are given in Appendix [A]. The result is that (for abelian gauge fields) take the 

form 

4 2) = ^((n„x)n^gF^ + x^[(x ± d)gF™])^Z 

+ \ * K 1^d c *M ^ + ^ ^d c i] >A T ^ (59) 



4? = l^TT \ m - D + %J P^ ~ — ~TT iJ P±* ) W °Q + f W-LcWc {(x ± D us ) q] 



+ "h.c". (60) 

The various applications of the equations of motion for £ are equivalent to the change of 
variables 

£ -> U + igx±A us - ^ {gx ± A ns ) 2 + l - (n_x) (gn + A us ) + ^ [ip d^A^A £. (61) 

The new collinear field £ on the right-hand side of (|6l|) can be shown to transform as 

£(x)^f/ us (x_)£(x) + O(A 3 0, (62) 

i.e. up to higher orders in A it transforms homogeneously with U us (x_) under the original 
transformation U us (x) that depends on all space-time components of x. The new effective 
Lagrangian is now manifestly invariant under the full set of gauge transformations, because 
ultrasoft gluon fields and space-time derivatives appear in covariant derivatives or field 
strength tensors only. Note that invariance of the Lagrangian under collinear and ultrasoft 
gauge transformations forces covariant derivatives in the x_ directions to always appear as 
in_D or in_Z) us but not as in_D c . 

The derivation of the corresponding results for non-abelian gauge fields is more com- 
plicated because of the interactions between ultrasoft and collinear gluons. This requires a 
careful treatment of the A expansion of the Yang-Mills Lagrangian which is left for future 
work. 



3.4 No renormalization 

The effective Lagrangian constructed so far reproduces the propagation and interaction 
of collinear and ultrasoft particles, but it does not yet include the effect of hard and soft 
fluctuations (loops) that contribute to a scattering process in full QCD. These effects are 
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included by matching perturbatively the QCD Lagrangian to the effective Lagrangian.[] In 
general, this modifies the coefficients of the various operators in the effective Lagrangian, 
and may induce new operators, so that 

A r ec = E^ — ' £ = £cvev- (63) 



We now show that the tree-level ultrasoft-collinear effective Lagrangian ( f4*8|) or (p^ ) is not 
renormalized at any order in the strong coupling, and at any order in the A-expansion, 
provided a factorization scheme can be used that does not break Lorentz invariance, such 
as dimensional regularization.Q In other words, no new operators are induced in higher 
orders in the matching and all coefficient functions retain their tree-level values, while the 
strong coupling evolves as in full QCD. 

Consider an iV-loop diagram in full QCD, where all loop momenta Zj are hard, which 
means that all components are of order A , and where all linear combinations of the Zj 
are hard. The external momenta (not necessarily on-shell) of the diagram consist of a 
number of collinear momenta pj, soft momenta Zv ms and ultrasoft momenta k nyUS . We 
include soft momenta here, because such hard diagrams can appear as subdiagrams in a 
larger diagram with soft loop momenta. Included in the definition of a hard subgraph 
is the instruction to Taylor-expand the integrand in all quantities that are subleading in 
A. All integrals are assumed to be defined in dimensional regularization. The resulting 
structure is a polynomial in external momenta except for the large components of collinear 
momenta. This precisely matches the structure of the operators that appear in the effective 
Lagrangian before integrating out also soft modes. The procedure outlined here defines 
a factorization scheme for the Wilson coefficients C,- and is similar to the scheme used to 



define non-relativistic effective theories in dimensional regularization [15 
The initial iV-loop diagram before Taylor-expansion takes the form 



N 1 

II dl i II 7^1 — ~7TTT -TTT ^ x polynomial, (64) 

j = X H ~r Zj Pj "T" 2-j K m,s + 2^ K n,us) 



where the product over k refers to the internal lines and the sums denote line-specific linear 
combinations of momenta with coefficients ±1,0, which we denote by L^^, Pk,c, Kk s and 
Kk,ns- Using the scaling rules for the various momenta, we find 

(L fc)h + P Kc + K Ks + K Kus ) 2 = L 2 kh + {n + P ktC ) (n_L fejh ) + O(A) . (65) 

The integrand of flBlj ) must be Taylor-expanded in all the terms of order A and higher, and 
we obtain a sum of terms of the form 

N l 

n " '< n m M+Mi „_ Li ^ x poly ' iomkl (66) 



4 We assume here that the strong coupling is small at the scales m and mX. When mX is of order of 
the strong interaction scale, only hard fluctuations can be integrated out perturbatively. 

It is assumed here without proof that dimensional regularization regularizes all hard and soft integrals. 
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with integers a k . But all these integrals vanish in dimensional regularization, because they 
could only depend on scalar products of the vectors [n + P kjC ) nt, which vanish because 
n 2 _ = O.f] It follows that the soft-collinear Lagrangian is not corrected by hard loops 
in the factorization scheme defined above. This unusual non-renormalization property is 
connected to the fact that, before fields are restricted to describe particular momentum 
modes, the collinear Lagrangian is in fact equivalent to full QCD as noted earlier, and 
acquires a Lorentz-invariant meaning only once it is coupled to a source that produces 
collinear particles. These sources are renormalized by hard fluctuations, but the Lagrangian 
is not. 

A similar reasoning can be used to show that soft fluctuations do not correct the soft- 
collinear effective Lagrangian either. Note that soft modes have virtualities of order A 2 , 
equal in magnitude to the virtuality of collinear modes. Nonetheless, soft modes can be 
integrated out, because they cannot appear as external modes under our basic kinematical 
assumption. A general soft diagram has all loop momenta Zj soft, which means that all 
components are of order A 1 , and all linear combinations of the Zj are also soft. The external 
momenta of the diagram consist of a number of collinear momenta pj and ultrasoft momenta 
k n . The matching onto the soft-collinear Lagrangian is again defined by Taylor-expanding 
the loop integrand in all small quantities and by evaluating all diagrams in dimensional 
regularization. With the same definitions as above the propagator denominators of a soft 
loop integral are 

(L k)S + P k<c + K k<us ) 2 = (n + P kfi )(n-L 



+ 



fc,s 

Lis + + PL + (n + P k , c )(n„K k>w ) + 0(A 3 ), (67) 



where the first term on the right-hand side is of order A and all others explicitly given are 
of order A 2 . For generic collinear momenta n + P ktC = 0(1), or else P m>c = if no collinear 
external momentum flows through line m. The soft matching correction then only involves 
integrals of the form 

N 1 1 

FT dd h TT 77 w s r- TT 7-5 s — x polynomial (68) 

with integers a k ^ m . These integrals are again zero in dimensional regularization. 

This argument does not cover exceptional configurations where collinear particles are 
nearly comoving, such that a subset of external collinear momenta sum up to nearly zero 
in an internal line n, i.e. n + P n)C ~ A. In this case the soft matching correction involves 
integrals of the form 

TV 1 

n^n 



i=i 



((n+Pn,c)(w-L n , B ) + L 2 )S + L ±n , s P ±n , c + P* c + ie) 



(i 



A more formal way to derive this goes as follows: perform a rescaling of all loop momenta, 



Under this scaling the loop integral transforms into itself times an irrational power of 9 for general d and 
integer a k - Hence the integral must vanish. 
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x TT — — - — r — — TT 7— — — x polynomial (69) 

with integers ak, m ,n, and these integrals are not obviously zero in dimensional regularization. 
We checked for some one- and two-loop diagrams that these integrals do vanish, because 
there is always at least one (n+Zj) integration where all poles lie on the same side of the 
real axis, but we do not have an all-order proof of this result. 

Under our kinematical assumption there are no tree diagrams involving soft particles. 
We therefore conclude that (barring the potential contribution from exceptional configu- 
rations discussed above) soft modes are not relevant for the construction of the effective 
Lagrangian. 

3.5 Heavy quark Lagrangian 

We will later discuss the production of collinear particles in the decay of a heavy quark 
near mass-shell. Heavy quarks near mass-shell are described by the Lagrangian of heavy 
quark effective theory (HQET), £hqet O [L7|, [18| . In our framework this Lagrangian 



corresponds to the interaction of heavy quarks with ultrasoft gluons. An interaction of a 
near-on-shell massive quark with a collinear gluon will throw the heavy quark off mass- 
shell and the corresponding process would have to be added to the effective theory as an 
explicit matching correction. The question we would like to address here is whether, when 
we combine the heavy quark with collinear and ultrasoft modes, the Lagrangian for the 
combined system is 

£ = £hqet + £scet, (70) 

where £scet is the soft-collinear Lagrangian, or whether the newly present collinear modes 
require that we add terms to the effective heavy-quark Lagrangian not present in the HQET 
Lagrangian with ultrasoft gluons only. Note that at this point we do not discuss sources 
(currents) which can convert a heavy quark into collinear particles. 

Before turning to this question we briefly review the standard HQET formalism for 
later use. Starting from the QCD Lagrangian with a heavy quark field Q, we define the 
four-component ultrasoft heavy-quark field Q v (x) = e vmv x Q(x), which varies only over 
distances x ~ 1/A 2 , m being the heavy quark mass. In terms of this new field the heavy 
quark Lagrangian reads 

£ = Q v (iIp-m(l-i))Q v = Q v (iip ns -m(l-i))Q v , (71) 

where to obtain the second equality we dropped the term containing the collinear gluon 
field, which vanishes by momentum conservation. In other words, the ultrasoft heavy 
quark does not interact with collinear gluons in the effective theory. One then performs 
the projection onto large and small components, 

fay ~ Qv 5 R-v ~ Qvi (72) 
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Figure 2: Collinear interaction of the heavy quark. 



and integrates out the small field H v . This gives the tree-level effective Lagrangian 



HQET 



h v 



ivD us + 



2m 



h 



and the leading-order equation of motion, 

ivD us h v = 0(\ 4 ) h v . 
Using ([74] ) Q v can be expressed as 

Q v = h v + H v = 



+ ^ + 0(A 4 ) h v . 
2m I 



(73) 



(74) 



(75) 



One also derives from the propagator of h v that h v ~ A 3 and H v ~ A 5 . Under ultrasoft 
gauge-transformations, h v — > U us h v , but h v — > h v under collinear gauge transformations. 

Now consider the interaction of a heavy quark with a collinear quark through two 
collinear gluons as in Figure where it is implied that both external heavy quark lines 
are near mass-shell. As a consequence of this the large components of the collinear gluon 
momenta must add to something of order A 2 . The intermediate heavy quark propagator is 
off-shell, and at first sight it might appear that this results in a new effective interaction of 
h v with collinear gluons or quarks which is not present in the HQET Lagrangian. However, 
diagrams of this type do in fact not represent permitted configurations for the physical 
process to which the effective theory applies. 

To see this, recall that the effective theory diagrams correspond exactly to those config- 
urations of momenta where the QCD diagram develops on-shell singularities. According to 
a theorem by Coleman and Norton these singularities correspond to classically allowed 
scattering processes. A heavy quark cannot emit or absorb two collinear particles without 
going off-shell, so that in such a classical scattering process one collinear particle must be 
in the initial and the other in the final state. In the effective theory for the decay of a B 
meson into collinear particles through a weak interaction current, there are no initial-state 
collinear particles. The invariant mass of the initial state is close to the mass of the heavy 
quark, and this excludes collinear particles in the initial state - B mesons consist only of 
ultrasoft heavy quarks and ultrasoft light quarks and gluons. We note that the same is true 
for soft particles: even if they appeared as degrees of freedom in the effective theory, there 
would be no vertex between a heavy-quark line and two or more soft quarks or gluons. 
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Figure 3: Tree-level matching of the heavy quark current. 



The conclusion of this discussion is that the combined heavy quark and soft-collinear 
system (in the absence of external sources, i.e. heavy-to-light currents) is described by the 
sum of the two Lagrangians ( |70|) with no additional interaction terms. 

4 Heavy- to- light current to order A 2 

In this section we discuss the representation of weak currents ipTQ, where T denotes a Dirac 
matrix, in the effective theory. We consider the kinematic region where large momentum 
is transferred from the heavy quark to the final state light quarks and gluons. 

4.1 Derivation of the effective current 

The purpose of this subsection is to derive the effective heavy-to-light current at tree level, 
including corrections of order A and A 2 . This extends previous work 0, where the effective 
current has been derived at order A . The emission of a collinear gluon from the heavy 
quark puts the heavy quark off mass-shell by an amount of order m 2 . The corresponding 
diagrams with off-shell heavy quark lines are not part of the effective theory, but must 
be reproduced by the effective current. Once a collinear gluon has been emitted from 
the heavy quark, the heavy quark stays off-shell when further collinear or ultrasoft gluons 
are radiated until the heavy quark decays into a light quark.f] The infinite number of 
tree diagrams that correspond to the emission of collinear and ultrasoft gluons has to be 
summed into the effective current as sketched in Figure |3|. 
In position space we obtain 

oo „ 

Jqcd(x) = $(x) r ^ / dzi ■ ■ ■ dz n D F (x - zi) ig4( z i) ' • ' D F (z n _i - z n ) igjf^(z n ) Q(z n ) 

n=0 

°° 1 1 

= fa) r £ ( - sm) ■ ■ ■ w ^ ( - jAW) q(x) 

7 At tree level, we do not have to consider the emission of soft gluons, since soft gluons do not appear 
as external lines for a physical process. 
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where Dp denotes the heavy-quark propagator and the covariant derivative — igA^ 
contains collinear and ultrasoft gluons, A = A c + A ns . The collinear gluon field next to 
the heavy quark field arises because the gluon at the beginning of the quark line must be 
collinear to put the quark off-shell. We omitted the ie prescription in the heavy quark 
propagators, since they are far off-shell by construction. The field Q v is defined in (fTSj) . 
The light-quark field will be discussed below. In the expression (|76|) gauge invariance of 
the current is not evident, since the collinear gluon field A% appears explicitly. Formally, 
gauge invariance can be made explicit by writing 

gA,Q v = [iip-m(l-i)]Q v -[iip ns -m(l-i)}Q v = [tip - m (1 - j)} Q v (77) 



where we used the equation of motion following from (T7T 
Our first task is to expand 



Q 



1 



Hp-m(l-if)) 



gA\^ I Q v 



(78) 



in powers of A. As before, we organize this expansion such that every term is covariant 
under the gauge transformations (0). For the discussion of reparameterization invariance 
it will be convenient not to use the canonical definition t> M = (n+ + n^L)/2 but to keep v as 
an independent vector. Consistency with our power counting scheme requires only that 



n+v ~ 1, ri-V ~ 1, v± ~ A 



(79) 



together with v 2 — 1, ny_ = 0, and n + ra_ = 2. To expand Q, we decompose the collinear 
gluon field into its different components and expand the inverse differential operator in 
ilp± and in^D. We do not (yet) separate the smaller ultrasoft gluon field in expressions 
such as (in + A) or A±, since this would not allow us to maintain invariance under the 
transformations (|Hp. The leading term in the inverse differential operator can be written 
as 

I 1 (l + i 4-\ 



Sn 



n_v 



in + D 2m J 



0) 



(in + D) ^_/2 - m(l - i) 

where here and in the following the inverse (in + D)~ 1 is understood to be defined with the 
+ie prescription. We then obtain 



Q 



Sq g n + A, 



ip±) S g n + A c — + gJji^ c 





n ~T) S gn + A c — + gn_A c — 
(-ifl±_) Sq g n + A c — + g4 Lc 



h v 

h v + 0{X 3 h, L 



(81) 
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Defining 

T = -^-gn + A c , (82) 
and using some Dirac algebra, the above expansion is rewritten as 

Q = (l-T )h v +(l- — (l+^T )f^h v 

\ ri-V 2 J 2m 

1 yi- / \ 1 fi- fi 



( (-ip±) T + gfl± c )h v -t ( (-m_D) T + g ri-A^)h % 

V / n. ii Am. v. \ / 



n_t> 2m ^ ' n_f 2m 2 

— (2< + (n+v) n») ( (-«£>„) T + <?,0 ^ 
n-V v 7 m + iJ v ' 

1 1 1 



({-i]) ± )T + gjL Xc )h v + 0{X 6 h v ). (83) 
n_f m m+D v 7 

This simplifies considerably when one uses ( p6|) to expand Tq^ in A, and then the identity 

(-i£P) (1 - W^ f ) + = iDWZ 1 - i£>£. (84) 

The term involving v± can be eliminated because 

1 / \ 2v 

[2v± + (n + v)n_) = n + . (85) 

n_v v ' n_v 

Combining various terms of similar structure and making use of the property (21) of WZ\ 
we find 



Q = wztfi + f^K- —^(iy>wzi-wziqpj)h v 

\ 2m J n_v 2m v ' 

- -L -L. ( Un,D + W# - 2ivD m - ft. 

n_v m + D \ l m J m J 

+ 0(X 3 h v ). (86) 

Neglecting corrections of order A 3 or smaller, the last two terms involving ultrasoft covariant 
derivatives can be modified by inserting WZ^ to the left and by replacing ip± by ip± us - 
This gives the final result 



Q = Wtf ( 1 + p±) h v - — t Upwtf - WZ^ ipj K 

\ 2m J ri-V 2m K ' 



1 1 t<2ivD + ^±\wZt-WZ< 2ivD m + *!-*!- 1 K 



in + D \\ m j [ m 

+ 0(A 3 /i v ). (87) 
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Note that the term with the derivative 2ivD = (n^v)(in + D) + 2iv±D± + (n + v)(iri-D) 
is of order A 2 , since the apparent order A term, in + DWZ\ cancels in the combination 
in + DWZ^ — WZ'in+D m = 0. Hence fl57|) does give a valid A expansion of Q. It is also 
manifestly collinear and ultrasoft gauge- covariant. In the absence of the collinear gluon 
field Q reduces to Q v as it should.^ 

Before we continue the discussion of ipTQ using (p7[) , we briefly mention an equivalent 
derivation of (|87D . Rather than summing tree diagrams we solve the Dirac equation 

(iip-m)Q = (88) 

in the presence of collinear gluon fields A%, together with the boundary condition that in 
the absence of a collinear field the solution reduces to the ultrasoft case (fT5|), 

Q{x) -> e- imvx Q v {x) forA£->0. (89) 

This equivalence can be seen from inserting (|77|) into ( [78|) , which implies that e~ imv ' x Q(x) 
is determined by the solution of the Dirac equation fl88|). We expand Q = + A + 
A 2 + . . . and solve the Dirac equation 



^— in + D - m (1 - if) J + + ^± f n _£) 



Q(o) + A Qd) + A 2 Q (2) + 1 = (9Q) 



iteratively in A subject to the boundary condition (|89|). The advantage of this procedure 
is that collinear and ultrasoft gauge invariance is manifest at every step of the derivation. 
For instance, to order A the equation is 



^ (m + D) - m (1 - ^ Q (0) = 0, 



(91) 



which is solved by = WZ^h v up to terms of order A 2 and also satisfies the boundary 
condition up to terms of this order. Substituting this into fl9"D| ) yields the equation for 
We solve this equation and satisfy the boundary condition by adding terms of the 
form WZ^ / us , where / us is an ultrasoft field. When this program is carried out to order A 2 
one reproduces the result (|8T| ) for Q. In this approach the particular form f(D^) WZ^ — 
WZ^ /(-Du S ) of many terms in ( |37D can be traced back to the requirement that Q — > Q v is 
demanded in the absence of collinear fields. 

Returning to the derivation of the effective current, we now have 

e imv - x Jqcd = J = (C + V + q) TQ (92) 



with Q given by (q7]). Inserting for r\ the expression that corresponds to ( p3|) and expanding 
this to the required order in A does not give the correct result. The reason for this is that the 



8 Although the last term in curly brackets in ( J87| ) is of order A 4 when A c is non-zero, it becomes of order 
A 2 as A c — > 0, because then (in+D)^ 1 — > (i7i + £> us ) _1 ~ 1/A 2 . However, this term exactly cancels against 
the first curly bracket when A c — > 0. 



24 



effective current is constructed such that it reproduces the on-shell matrix elements with 
a current insertion in QCD. In the effective theory this amounts to calculating insertions 
of the effective current together with time-ordered products with terms of higher order 
in A in the effective Lagrangian. Different forms of the effective Lagrangian which are 
equivalent by the equations of motion give different results for the effective current beyond 
leading order in A. The correct current is automatically determined by matching the matrix 
elements in the full and in the effective theory. In order to avoid the explicit calculation 
of the time-ordered product terms, we couple the current to a weak external field B and 
add the terms 

Cj = e~ im v ' x $T QB + "h.c." = (£ + fj + q)j + "h.c." (93) 
to the Lagrangian (|3"2"D. We now repeat the derivation of the effective Lagrangian (ESI) in 



Section 3.2 and follow the modifications induced by the presence of the additional source 



term. The solution for the field rj in (^) now reads 

V = — ^ ^ + gM + j)- (94) 

Inserting this into (|32|) and ( fffi) results in 

C 3 = (V + q - (e ilp ± + qg4c) (infi] 

= (V - f 0± (*n£) ~ 1 ^ + q~q ^ (1 - ZW^ j + "h.c." (95) 




where we expanded in A up to second order, used (|36|) and neglected terms quadratic in 
the weak field B. In Section we used the £ equation of motion to pass from ([TI]) to 



With the term Cj added, the leading-order equation of motion ([42]) for the collinear 
quark field £ is modified to 



so that in passing to (|43|) one picks up an additional contribution 



q(l -ZW ] )^^j + "h.c." (97) 



to Cj. We then obtain the result 



Cj = (t - e $ ± (m + "D) ~ 17 y + qZWU J + "h.c." 



(9£ 



for the source Lagrangian in the effective theory with the Lagrangian (fJ8j), which implies 
the effective current 



J = h-£ L (m+&) 1 4± + qZW^j TQ, 



(99) 



25 



where it is understood that derivatives and their inverse do not act on the overall phase 
factor e~ imvx . The field redefinition 



£ - £ - (1 - Wtf) 



(100) 



implicit in using the field equation for £ in (|4l| ) has rendered this expression for J mani- 
festly invariant under the gauge transformations (|HD , contrary to the original expression 
(|92|). Under a gauge transformation the external current transforms as j —>■ Uj, which is 
compensated by the transformation of the round bracket in fl98|). This is similar to the 
structure of the effective Lagrangian, where ((48|) is invariant but the starting expression 
(|32|) is not. In other words, it is the field £ after the (implicit) field redefinition ( |100| ) that 
transforms according to the gauge transformations (|14]) but not the original field £. 

The ultrasoft quark field does in fact not contribute to the heavy-to-light current at 
relative order A 2 . Since q scales as A 3 , we have 



qZW ] TQ « qZW^T WZ ] h v = qT h v ~ 



(101) 



for colour-singlet currents. The final term qTh v can be set to zero, because we assume that 
there is large momentum transfer to the final state, so that the operator must contain at 
least one collinear field in order to contribute to such final states. Inserting the result for 
Q from ( |37D into (|99|), we obtain for the current including corrections of order A 2 : 



$(x)TQ(x) 



QCD 



J(0) + j(D + J(2) + 



(102) 



with 



jC 1 ) = -£0 . (in.fi) ll ^TWZ ] h v - —tT^- WiiPsWZ^] h v , 
v ' 2 n_v 2m LL JJ 



(103) 



= lTWZ^h v - 

2m n_v 



1 .^r^ - ^ 



ri-V 



Am 



in. 



,DWZ ] 



ri-V in+D 



2ivD + 



m 



where the double-bracket is defined by 

~f(D)WZ^}} = f(D)WZ ] - WZ^f{D v 



(104) 



The double-bracket structures generically appear in (|103|) , because they guarantee that the 
power counting of the various terms is preserved when the collinear gluon field is taken to 
vanish. Each double-bracket vanishes separately in this limit, so that terms with l/(in+D) 
are never enhanced. 
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4.2 Reparameterization constraints 

Contrary to the soft-collinear effective Lagrangian, the effective heavy-to-light current does 
receive radiative corrections, which may induce non-trivial Wilson coefficients multiplying 
the operators that appear at tree level, as well as new operators with Wilson coefficients 
whose expansion begins at higher orders in a s . We discuss in the following how the various 
operators are related by reparameterizations of the velocity vector v and the light-like 
vectors n±, which can be used to restrict the number of independent Wilson coefficients. 
Such reparameterizations of the soft-collinear effective theory have been used in EU 



21| to constrain the renormalization of some of the power-suppressed operators in the 
effective Lagrangian. In Section |3.4| we have seen that none of the possible operators in 
the effective Lagrangian can acquire a Wilson coefficient different from its tree-level value, 
so reparameterization invariance does not provide further information on the Lagrangian. 
In the following we show how reparameterizations relate the various terms of the current 
( |I03|) . We do not attempt a complete analysis, which would have to treat the set of all 
possible operators at order A and A 2 , not only those that appear at tree- level. 

The effective current depends on the arbitrary velocity vector v defining the projection 
h v of the heavy quark field. As in heavy quark effective theory the current must be invariant 
under a velocity transformation v — > v' = v + Sv with (V) 2 = I, where 8v ~ A 2 , since 5v is 
related to the ambiguity in splitting a heavy quark momentum in its kinematical piece mv 
and a residual momentum of order A 2 . Invariance of the QCD current implies that J(x) 
or, equivalently, Q(x) should transform as 

Q(x) -> e imSv - x Q(x). (105) 



Working out the transformations of the various terms in (p7|) , we find that velocity repa- 
rameterization invariance ties 

Wtf^K, (106) 
2m 

which appears at order A 2 , to the leading-order term WZ^ h v just as in heavy quark effective 
theory, but no other relations are enforced at order A 2 . Invariance under velocity reparam- 
eterizations implies, however, that the two terms in the definition of the double-bracket 
( |I04j) should not be separated. 

Manifest reparameterization invariance is achieved by expressing all operators in terms 
of the field Q v defined in (|75|) and the differential operator 

yM = v» + (107) 

m 

rather than h v , and iD 11 |22J, since 

Q v (x) - e m5v - x Q v (x), /(V) Q v (x) - e imS ™ /(V) Q v (x). (108) 
We can indeed cast the expression ( |37D for Q into this form, the result being 

n-V 2 11 JJ 
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V 2 - 1 



{n_V n + V + f ± f ± } WZ^]]Q v + 0(\ 3 Q V ). (109) 



Here we have used that up to higher orders in A we can approximate l/(n_v) ~ l/(n_V) 
and m/(n_vin + D) « 1/(V 2 — 1). This can be written even more concisely by noting that 

nJVn + V + Y±y± = yy-in+D^^n-V + yA +C(A 3 ) (110) 
under the double-bracket, so that 

q = wz^q v - -J— \\yywz^]]Q v + o(\ 3 q v ). (in) 

By the same line of reasoning the light-like vectors n + and n_ that characterize collinear 
particles in the effective theory are not uniquely determined. The allowed reparameteri- 
zations that are consistent with the power-counting scheme for fields and derivatives, and 
leave n + n^ = 2 and n 2 ± = invariant, can be divided into three classes that correspond to 
particular Lorentz transformations of the vectors n + and n_ 0, 11, 14]: 



a) longitudinal boosts: n + an + , n_ — > a -1 n_, 

b) "minus-transverse" boosts: n + — ► n+, n_ — > n_ + 2ej_ — n + , (112) 

c) "plus-transverse" boosts: n_ — > n_, n + — >■ n + + 2/j_ — f]_n_. 

The velocity vector is considered as independent of n± and remains fixed under all three 
transformations. Here e± and f± are arbitrary vectors in the transverse plane, and power- 
counting for n_D, D± and v±_ requires a ~ 1, and e±, f± ~ A or smaller. Of course, each 
single operator in the effective theory must still be manifestly invariant under rotations 
within the transverse plane (which leave n + and n_ unaffected). 

The requirement that a vector a remains invariant under the above transformations 
allows us to deduce the transformation of its transverse, plus- and minus-components. 
Under longitudinal boosts a) 

n + a^an + a, n^a —> oT l n_a, a± — > a_i_. (H3) 

Under minus-transverse boosts b) 

n + a — > n + a, 

n_a — > 7i-a + 2ej_a — e\ n + a, 
a ± ~^ a 1 ~~ n + ae i — e±an+ + e 2 ± n + an+, (H4) 

while under plus-transverse boosts c) 

n + a — > n + a + 2f±a — f]_ n_a, 
n_a — > n_a, 

— > — ri-a f±_ — f±ant + f]_n + an ! t. (H5) 
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The transformation of £ follows from the invariance of the full spinor ifj c and the definition 

£ = (Ji_ifi + ) ipc/^i which gives 



116) 



The Wilson line WZ^ is invariant under the transformations a) and b) and transforms 
under plus-transverse boosts as 



wz ] -> wzt 



in+D 



if±D ± WZ j }} +C(A 3 ). 



(117) 



We now discuss the implications of the different reparametrizations of collinear fields in 
turn. Longitudinal boosts constrain the allowed operators in the effective theory, but do 
not relate different operators. It is easy to see that for each operator in the effective heavy- 
to-light current ( |103| ) the powers of a induced by the rescaling of n + and n_ compensate 
each other, so each term is invariant. 

In contrast, transverse boosts lead to non-trivial relations between different operators in 
the heavy-to-light current. In particular, the leading order current £ V WZ^h v is not invari- 
ant under transverse boosts. It is easy to see that the transformation of £ is compensated 
by that of 

which suggests that we express the current in terms of 



1181 



(119) 



which is simply £ + fj in the absence of the ultrasoft quark field. This relates the first term 
of J*- 1 -* in ( |103| ) to the leading-order current and the second term in to the second term 
in JW. To exploit transverse boost invariance further, we write 



as 



J = ETWZ^Q V - — Hrp 

Ti-V 2m 



in*D 



2ivD WZ ] 



ip ± wtf]]Q v - J_sr%^±" 

ri-V 4m L 
1 WiJPxiJP 



in+D 



in_DWZ^ 



rn 



Qv 

Q v + ... 
(120) 



where the J, are defined by the successive five terms in the expression for J. 

We consider first the plus-transverse boosts c). The leading-order term Jo is not invari- 
ant due to the transformation ( |117| ) of WZ*. The non-invariant term is cancelled by the 
transformation of the term J3 containing the derivative 2ivD, so the Wilson coefficients 
of these two operators must be related. We note that if we had written 2ivD in terms 
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of in + D, iri-D and 2iv±D } reparameterization invariance would have tied all these terms 
together into 2ivD. Plus-transverse boosts yield no further relations up to order A 2 , since 
the variations of Ji,2,4 are all of order A 3 . 

The leading term Jo is invariant under minus-transverse boosts b). The order- A term 
J i is not invariant. Its variation reads 

8 h) J x = — ST — [[{2 tf ± ip ± - fi_fi + ie ± D} WZ*]] Q v + .... (121) 

Ti-V 2m 11 JJ 

The first term in curly brackets is cancelled by the transformation of J4, the second by the 
transformation of J2. This implies a relation among the Wilson coefficients of these terms. 
We note that if we had decomposed ip±ip± into its symmetric and anti-symmetric part, 
both parts would have taken part in the cancellation of 5^ J\. 

We conclude from this that for a given Dirac structure T the operators in the effective 
heavy-to-light current that appear at tree level are divided into two sets that correspond 
to { J , J 3 } and { Ji, J 2 , J4} in ( |120| ). The Wilson coefficients of the operators in each set 



are related by reparameterization invariance. The fact that the operator structure J\ is 
not related to J under transverse boosts has already been noted in [^3[ . 

Upon renormalization we may also encounter new power-suppressed operators, which 
are not present at tree level. These operators will also be related by reparameterization 
symmetries, but the derivation of the complete set of relations is beyond the scope of the 
present paper. Note that the terms J3 and J4 involving 2ivD and ilp±ilp_\_/m, respectively, 
are not separately invariant under heavy quark velocity reparameterizations, but their 
sum is. This suggests that there may be a further relation that connects the two sets of 
operators defined above, but since the non-invariant term in J 3j4 is of order A 3 , one would 
need to analyze the set of higher-order operators to draw a definite conclusion.^ Having 
exploited the reparameterization invariances we will from here on present our results with 
the particular choice v = (n + + nJ)/2 for the heavy quark velocity vector. 



4.3 Expansion in A 

The expansion of the effective heavy-to-light current in terms of operators with definite 
(homogeneous) A scaling is performed in the same manner as for the collinear Lagrangian 
in Section |3.3| . The derivatives and inverse derivatives in the current that operate on £ to 
the left are expanded as 




-_l_us 




'122) 



where A us is now evaluated at x = x_ . Similarly, the quantity Q in fl57|) is expanded using 
(p5p for WZ\ and all ultrasoft fields are multipole-expanded according to ([5ID. Multiplying 

9 This relation becomes even more suggestive when the current is represented as STQ with Q expressed 



by (111). Then the only quantity not manifestly invariant under transverse boosts is WZ^ , and the 



variations of the two terms in (111) compensate each other 
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everything together we obtain the equivalent of ( |103| ), where now every term has a definite 
scaling in A. To present the result we find it convenient to split the heavy-to-light current in 
the effective theory into two terms (denoted and J^) which are separately invariant 
under the reparameterizations of n± discussed in the previous subsection. Neglecting terms 
beyond order A 2 , we obtain 



'■4>{x)TQ{x) 



QCD 



-«{j(A) +J (B)j 

{ J (A0) + J (A1) + J (A2) + J (m) + J {B2) } (123) 



with 



j(A0) 
J(A1) 

j(A2) 



lTW c h Vl 

lYW c [x ± dh v }-tity Xc (in + "D c )~ lr t±TW c h v , 
er W c i\{n-x) [(n+d) h v ] + X - Xfl± x u± [d»d v h v ] + 



2m 



in.D r 



in_DW c h v + £T 



W c , gn + A u 



in + D c 

- £iip ±c (in + %) ~ l7 t±TW c [x ± d h v ) + I g4 ±ua (in+%) ~ l7 t±TW c K 



(124) 



and 



J(B1) 
j(B2) 



-iv-^iip Xc w c h v , 

2m 

7ji_ . - d. 



(125) 



-I r Xfl± ilp ±c W c [& 1 h v }-lvj-{ [ g A ±us , w t 



2m 



+ 



in-D W r 



1 Uppity. 



in + D c 



m 



W c h v + tilj) ±c (in+% 



2 2m 



where on the right-hand sides the ultrasoft fields h v and A ns are evaluated at x = x_, 
and derivatives on ultrasoft fields only act within square brackets. The leading-order 
current J^ ) coincides with the result derived in by summing momentum space Feynman 
diagrams. When we compare our terms at order A with the result obtained in p0[ , we find 
agreement for the term j( A1 ^ (except for the first term coming from the multipole expansion 
of the heavy quark field), but find that the term J( B1 ^ has been missed there. The second- 
order corrections in j( A2 ^ and j( B2 ^ have not been considered before. 

The A expansion of the heavy-to-light current as given by ( |124j ) and ( |125| ) corresponds 
to the effective Lagrangian fl5"2"|). Gauge invariance is therefore manifest only for ultrasoft 
gauge transformations restricted to depend only on x_. The same steps that led us (for 
abelian gauge fields) from the expanded Lagrangian ( |52"D to the Lagrangian with terms 
(j58|)-(l60D, automatically convert the current into an expression where ultrasoft gauge fields 
appear only in covariant derivatives and field-strength tensors as in the corresponding 
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Lagrangian. This is derived most conveniently by noting that the manipulations of the 
Lagrangian correspond to the field redefinition of the collinear quark field £ given in (|6l|). 
Inserting this into the source term for the heavy-to-light current in (|98|) we find that partial 
derivatives and ultrasoft gauge fields in the heavy-to-light current given by ( |124j ) and ( |1 25|) 
recombine into covariant derivatives. The currents J^ A0 ^ and J^ 1 ^ in ( |124j ) and ( |125 ) 
already assume their final form. For the other terms we obtain (for abelian gauge fields 
only) 



J(A1) 
j(A2) 



jm 



iTW c [x ± D w h v ]-&l? 



±c 



in+*D c ) 1 ^±FW c h v , 



I 1 

[W h v ] + 



2m 



in+D r 



in_D W r 



4- 



m. 



±c 



in + D c 



m 



W c h v + iiip 



in + < D c 



(126) 



2 2m 



This form of the current corresponds to the manifestly gauge-invariant effective Lagrangian 
given by fl58|)-(|5"0D for abelian gauge fields. The non-abelian case is left for future work. 



5 Form factors for heavy-to-light meson decays 

We now apply soft-collinear effective theory to the heavy quark expansion of transition form 
factors for B decays into light pseudoscalar or vector mesons in the kinematic region where 
the energy of the final state meson is of order of the heavy quark mass m. Because of the 
large momentum transfer to light partons, collinear modes are relevant in this kinematic 
region. 

The transition form factors encode the hadronic effects relevant to semi-leptonic decays, 
but also provide important non-perturbative input to rare and non-leptonic B decays in 
the QCD factorization approach f24] , |25| , |26| . In QCD one defines three independent form 
factors for B decays into pseudoscalars, and seven independent form factors for decays 
into vector mesons.^] At leading order in an expansion in 1/m, the number of independent 
form factors for B decays into light pseudoscalar (vector) mesons is reduced to one (two) 
0, |2"7fl . This leads to form factor relations which reduce the hadronic uncertainties in the 
calculation of exclusive B decays. 

In H two of us have shown (to first order in the strong coupling a s ) that the form 
factors factorize at leading order in 1/m into a soft-overlap contribution multiplied by a 
short- distance coefficient Cj and an additive hard spectator-scattering correction A/j. For 

10 We use the same definitions as in [Bl. They are summarized in Appendix H for convenience. 
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instance, the three form factors fi{q 2 ) describing the decay into a light pseudoscalar meson 
can be written as 

fM 2 ) = C?\E)^\E) + ^Af i (E), (127) 

where the energy E of the final state meson is related to the momentum transfer q 2 by E = 
(m 2 B + m 2 P — q 2 )/{2rriB) and £p is the universal "soft" form factor defined in the effective 
theory which parameterizes the soft-overlap contribution at leading order in 1/m. The tree- 
level coefficients Cf^ are given by calculable kinematical factors. Radiative corrections to 
and the hard scattering corrections A/j have been obtained in in a "physical" 
factorization scheme defining the soft-overlap form factors. The radiative corrections to 
Cf^ have also been computed in soft-collinear effective theory with MS subtractions. 
The implications of the two calculations coincide for QCD form factor ratios from which the 
soft-overlap form factors drop out. The effective theory approach provides the appropriate 
tools to prove a factorization formula such as (|127|) to all orders in a s . A first step has been 
taken in H, where it was shown that the relations among the soft-overlap contributions 
to the QCD form factors hold to all orders in the absence of hard spectator interactions. 
To complete the proof one would have to demonstrate that the infrared contributions that 
appear in the calculation of the hard scattering term can be absorbed into the universal 
soft-overlap form factor to all orders in a s . 

Combining the soft-overlap and hard scattering terms requires an extension of the ef- 
fective theory that has not been formulated so far. In the following discussion we focus on 
the soft-overlap contributions beyond leading order in the 1/m expansion. These are con- 
tributions where the spectator quark in the B meson does not undergo a large-momentum 
transfer interaction, so that the final state consists of the ultrasoft spectator quark together 
with other small momentum modes, and a cluster of partons with large momentum that 
evolves out of the b —>■ light quark transition. The invariant mass of this cluster of partons 
is of order mA, where A is the strong interaction scale, since the interaction of an ultrasoft 
gluon with momentum A and a nearly on-shell quark or gluon with momentum of order 
m puts the large momentum parton off-shell by an amount mA{£] This implies that we 
should take the scaling A ~ (A/m) 1 / 2 for the expansion parameter in soft-collinear effective 
theory. We may then write 

h{q 2 ) | soft overlap = Cf £<°> + £ + Cf *f + ...), (128) 

3 

where the functions Q are defined by matrix elements of the terms of order A and A 2 , 
respectively, in the effective current derived in Section [|. The purpose of this section is then 
to see whether form factor relations may survive when one includes power corrections. Note 
that we have to go to second order in A to cover 1/m corrections because A ~ (A/m) 1 / 2 . 

Before we begin the discussion of power corrections it is instructive to see how the 
multiplicative structure Ci(E) ^i(E) in the factorization formula ( |127| ) arises in the effective 

11 From this point of view the configurations which are considered in the hard scattering approach are 
special, since they consist of only large momentum partons, in which case their invariant mass must be of 
order A 2 , the light meson mass squared. 
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theory formalism. With hard radiative corrections included, the leading-order matching 
equation 

'^(x)TQ(x) 

is modified to 

^(x)TQ(x)] QcD = e- mv - x J2 f ds& T (s, f j,)O i (x + sn+). (130) 

Here we defined a more general set of current operators in the effective theory 

Oi(x) = Z(x)TiW c (x)h v (x-) , (131) 

where Fj are appropriate Dirac matrices, which may be induced by radiative corrections, 
and Cf (s) are the corresponding hard coefficient functions. The integral over s appears, 
because the component n + p of collinear momenta is of the same order as the hard momenta, 
so that renormalization in coordinate space is non-local.Q Choosing x = we obtain for 
the matrix elements that define the heavy-to-light form factors 

(M(p')\$TQ} QCB \B(p)} = W° daCi,{s,^(M(jpf)\e^^r i W e K]e-^ p \B(v)) 

i J -°° 

= Y^C l r (n +P '-(m B -m)-ie,fi) (M(p')\ [^W.K] \B(v)), (132) 

i 

where we introduced the translation operator e mP and used P tM \B{v)) = (mB — m)v^\B(v)) 
on the heavy meson state in the effective theory. We may replace n + p' -(mg-m) by n + p' 
in ( |132|) , since the difference is a higher-order effect in the 1/m expansion. The momentum 
space coefficient functions C r (n + p') defined by 

C l r (n + p',fx) = f° dse isn+p ' ' C l T {s^) (133) 

correspond to the coefficient functions computed in 0. The matrix elements on the two 
sides of ( |132j ) define the form factors in QCD and in the effective theory, respectively, so 
that ( |132| ) displays the multiplicative structure of the matching equation for the soft-overlap 
form factors. 

12 The relation between the renormalized and bare operators in the effective theory is also of this non-local 
form, 

Oi(x,a,ii) = J ds' Z(s- s',fi) 0^ rc (x,s'). 

The non-local contributions to the Z-factor come from the ultraviolet divergences of collinear loops. Ul- 
trasoft loops give local contributions proportional to S(s — s'). 



: Z(x)TW c (x) h v {xJ) 



(129) 
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5.1 Leading power 



We first recapitulate the derivation of the form factors at leading order in A, where the 
effective current is given by £ TW c h v according to fljjjj ). Because of the projection prop- 
erties of the light collinear field, 4-£( x ) = 0, and the heavy quark field, (1 — i/i) h v (x) = 0, 
the most general decomposition of the matrix element of the leading-order current between 
a B meson and a light pseudoscalar (P) or vector meson (V) reads 



(P(p')\trW c h v \B(p)) = tr 
(V(p',e)\trW c h v \B(p)) = tr 



A (o) 4+4- r 1 + i 



A 



(0) 



4+4- th i + 



(134) 



where the current operator is evaluated at x — 0. In the following we consider only 
pseudoscalar B~ or B° mesons. Note that the dynamics of the fields £ and h v in the 
effective theory, defined by leading order parts of the Lagrangians fl3"T| ) and ([73]), does not 
change the projection properties of these fields. Without loss of generality we choose the 
four- vectors v, n + , ri- to satisfy v = (n + + nS)j2 as before, and in addition p = uibv 
and p' = \{n + p')n_ + |(n_p')n + . Up to relative corrections of order m 2 PV /E 2 ~ A 4 
we have p' = En_, where mpy is the mass of the light meson. The universal functions 
Ap\E, v, n_) and Ay (E, v, e) can be decomposed into their most general independent 
Dirac structures allowed by Lorentz invariance and parity, 



aP 

|(0) 



(E, v, n. 



2E^\E), 



-IE f ± 75 d° } (E) -2E(ve*) 75 (E) 



(o). 



(135) 



leaving only one universal form factor ^ (E) for decays into light pseudoscalar mesons, and 
two form factors ^(E), £,| (E) for decays into transversely and longitudinally polarized 
light vector mesons. Here e± = e — |(n + e)n_ — |(n_e)n + is the polarization vector for 
a transverse vector meson. (We slightly changed notation compared to |§ and included 
the dependence on the polarization vector and factors 75 into the functions Ap and Ay.) 
The resulting form factor relations can be found in [S, |27|, or can be read off the leading 



contributions (involving only fp±l defined below) to the formulae given in Appendix |B[ 



5.2 Beyond leading power 

Beyond leading order in A one could continue to write down the most general decompo- 
sition of the matrix element of every term in the effective current (|103|) , as well as the 
decomposition of time-ordered products of the current with higher-order terms in the ef- 
fective Lagrangian. This would give rise to a large number of non-perturbative functions 
Q ' defined in the effective theory, and an expression such as ( |128| ) for the QCD form 
factors in terms of these functions. This approach would be useful if the form factors in 
the effective theory could be computed in a simpler way than the form factors in QCD, 
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for instance in lattice QCD. This is unlikely for soft-collinear effective theory, where the 
operators are all non-local in light-like directions. One must then construct combinations 
of form factors in full QCD from which the unknown non-perturbative functions drop out 
to a certain order in A. For the derivation of such form factor relations a simpler approach 
is available that avoids the decomposition of every operator and time-ordered product. 
We start by writing the matrix element of the heavy-to-light current in full QCD as 



(p(p')\i>rQ\B(p)) 

+ tr 



tr 



A (++) 4+4- r l_ 



4 



+ tr 



A 



(+-) 4+4- 



( - +) 4-4+ i 
a p — ^— i - 



+ tr 



A 



{—) 4-4+ 



4 

i -j 



(136) 



with an analogous definition for decays into vector mesons. Note that the decomposition 
of the right-hand side of ( |136|) simply uses the projectors (1 ± and rjbJj/i^/A on the 
large and small components of the quark fields, labelled by + and — , respectively, and is 



\{kl) 



completely general. The functions A p 

as 



and Ay l> with k,l = +, — can again be decomposed 



.4 



(hi). 



p (E,v,n. 



2Ef£ l) (E), 



A^(E,v,n.,e) 



2Ef ± 75 /f \E) -2E{ve*) 75 f\ m \E). 



{hi), 



(137) 



(kl) ( kl) 

Among the 4 + 8 form factors fp and f]_ y only 3 + 7 are independent. This follows from 
the equations of motion for light and heavy quarks in QCD and translational invariance, 



q»{P\j>^b\B) 
q»(V\j> W b\B) 



(m-m q ) (P\i)b\B), 
-(m + m q ) {V\$ lh b\B) 



(138) 



with q = p — p'. This implies 



m B -m (++) q , ( _ +) ( _ 
J P II T l J P II JP 



m B 



')-{*% 



+-) 



A-+) 



f P J ] } = 0(X 3 f^), (139) 



where we have neglected light quark masses m q and anticipated that f P + ^ is the leading 
form factor in the A expansion. Note that itlb — m ~ A 2 , mpy ~ A 2 , and that there is no 
relation for decays into transversely polarized vector mesons. The remaining independent 
3 + 7 form factors are in one-to-one correspondence to the conventional form factors pa- 
rameterizing heavy-to-light transitions for vector, axial vector, and tensor currents. Their 
definitions and the transformation between the two sets of form factors are given in Ap- 
pendix |B[ 

To find the form factor relations we take the effective current (|103|) and determine the 
order in A where one obtains a contribution to the form factors // (with i = P, 4_, ||) 
by decomposing T in each term into its four projections analogous to (|136|). For instance, 
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the leading-order current only contributes to A^ + ^ (with M = P,V) as seen in the 
previous subsection, so 

ft + \E)=^\E)(l + 0(a s ,X)) (140) 



for the leading form factor in the large-energy /heavy quark mass limit. The order A terms 
in the current contribute only to A^ + \ A^~' and A^^ but not to A^~\ Finally, the 
A 2 suppressed terms in contribute to all form factors. This can be summarized in the 
relations 

ff-\E) ~ Ajf +) (£), ft + \E) ~\ft + \E), ft~\E) ~ X 2 fl ++ \E). (141) 

In each case these relations specify the largest term in the A expansion. This tells us 
that there are no form factor relations when power corrections of order A 2 are included. 
However, at order A we can neglect // (E), which leaves us with 3+6 form factors for 
pseudoscalar and vector mesons, respectively. The relations ( |139| ) approximated to order 
A give two constraints, which implies 2+5 independent form factors. Therefore, we have 
1+2 form factor relations at order A. 

5.3 Form factor relations 

The three form factor relations that continue to hold including the leading power corrections 
can be found from the relations between ff l {E) and the conventional heavy-to-light form 
factors in Appendix [B], neglecting j i [E). The relations read 

/+ " h Q 2 



Rp = i±ZJl = g_ (l + Ofo.A 1 

f T m B [m B + m P ) v 




R± = 7 2~^\ = 7 — — , ( 1 + 0(a s , A 

/ q \ tub [mB + my) v 



V 



Rw 




t,-(i-4)t, 

V mi I 



= -^(l + 0(a s ,\ 2 )) (142) 

for form factors that describe the decays into light pseudoscalar, transversely and longitu- 
dinally polarized mesons, respectively. Note that the vanishing of the right-hand sides of 
( |142|) with the momentum transfer q 2 follows from the absence of massless poles in hadronic 



37 



matrix elements for q 2 —* 0, and is thus exact in QCD. Other form factor relations valid 
at order A do receive A suppressed power corrections in the effective theory, for instance 




V- (l + ^-Y A 1 = 0(\), (i-^ Tl -T 2 = ^0(A), (143) 
V m B J V m B m B 



where the first relation is proportional to f[_ + f± ^ and the second to f± + ^—f± ^ As 
observed in [28fl, the form factor combinations (|143|) describe the transition of a collinear 



quark with negative helicity into a vector meson with positive helicity. We explicitly see 
that helicity retention from quark to meson does not hold beyond leading order in the A 
expansion. In this respect our conclusion differs from the one in P0 |, where it is stated 



that the first equation in ( |143j ), which relates the vector and one axial- vector form factor 
for B — > V decays, does not receive order A corrections. (The tensor form factors fy, Ti 2 ,3 
were not discussed in [pp . 

If we include order A 2 corrections, non-trivial form factor relations do not survive. 
However, ([142) may still be useful if we restrict our analysis to small momentum transfer 



and impose the scaling q 2 ~ A 2 . In this case the undetermined A 2 corrections are always 
multiplied by q 2 , and can be neglected to the considered order. The same holds for the 
radiative corrections to ( |142| ). 

To see how the form factor relations are satisfied in a specific non-perturbative model, 
we compare the left-hand sides of (|142|) with the form factors computed in the QCD 
sum rule approach in 1211 RO] to the right-hand sides. To make this comparison more 



accurate we also include on the right-hand sides the known a s -corrections to the heavy 
quark limit, which can be obtained from Section 3.3 of ||. Interestingly, the hard-scattering 
contribution drops out from all three ratios Rp t ±,\\, and the perturbative correction is due 
to the coefficients in ( |127| ) alone. The ratio R± receives no radiative corrections 



at all at leading power in A because of the helicity argument given above. On the left- 
hand side of Fig. |] we plot the difference between the functions Rp,j_,\\ and their values in 
the symmetry limit (defined as the limit where power corrections and a s corrections are 
neglected). The dashed curve gives the QCD sum rule result, and the solid curve gives 
the result in the effective theory at leading order in A, including radiative corrections. 
The difference between the two curves is therefore a measure of A 2 power corrections, or 
corrections of order ct s \. Because of the q 2 suppression of the right-hand sides of ( |142| ) we 



expect only small corrections. Indeed, in all three cases the deviations from zero are at 
most 3% for q 2 up to 7 GeV 2 . 

On the right-hand side of Fig. [| we plot the ratio of the functions Rp,±,\\ and their 
values in the symmetry limit minus 1. In this case the q 2 suppression is divided out, and 
we expect deviations from of the order of A 2 and/or a s . In the case of Rp and R\\ 
the effect of radiative corrections is typically 10%. If we take the difference between the 
light-cone sum rule predictions and as a measure for the power-corrections in ( |142| ), we 
deduce effects of order 15-20% which is compatible with A 2 = AQ CD /m or a s \. Of course, 
at this point it is unclear whether the difference between the two sets of curves is a real 
power correction or an artefact of the QCD sum rule calculation. 
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Figure 4: Illustration of the relations (|142|) for B — > ic and B — >• p form factors. On the 
left-hand side we plot the ratios Rp,±,\\ minus the corresponding values in the symmetry 
limit. The right-hand side shows Rp,±,\\ divided by their symmetry limits minus 1. The 
dashed line corresponds to the form factors predicted from QCD light-cone sum rules in 
|29| , ^0| . The solid line is the order A result supplemented by a s corrections as calculated 
in ||. The tensor form factors are evaluated at the scale m. The grey error band reflects 
the theoretical uncertainty from varying the scale of a s from m/2 to 2m. 
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A subtle question concerns the power-counting of hadronic matrix elements in soft- 
collinear effective theory. In contrast to the well-known Isgur-Wise form factor in HQET 
31], which is constrained to scale as A by its normalization at zero recoil, the power 
counting for the universal form factors ^°\e) is not determined by the effective theory. 
The scaling behaviour of the quark fields and of the meson states with A is not sufficient 
to determine the corresponding scaling behaviour of form factors. Formally, we have from 
the normalization of hadronic states 

(PGt/)|~A- 2 , \B{m B v)) ~ A~ 3 , (144) 

where we used p' ~ En_. Together with £ ~ A and h v ~ A 3 this results in 

(P(p')\ZTW c h v \B(m B v)} ~ A- X X(A). (145) 

Here we introduced the function X(X) which parameterizes the probability to find the 
parton configuration corresponding to the current operator in the incoming and outgoing 
meson states. While the probability to find a nearly on-shell 6-quark in a B meson is 
of order one, the configuration with one collinear quark (jet) carrying most of the light 
meson's energy is expected to be suppressed. The suppression factor is usually estimated 
from the end-point behaviour of the meson's two-particle light cone distribution amplitude 
(see |§ and references therein), which gives 

{t§\E) , if{E) ,(ve) ^\E)} oc A 3 . (146) 

From this one concludes that X(X) ~ A 4 . In our discussion of power corrections to the tran- 
sition form factors we implicitly assumed that the matrix elements of higher-order terms 
in the effective current have the same A 4 suppression as the leading-order matrix element. 
There might be further dynamical suppression of the higher-order matrix elements, but an 
understanding of this requires knowledge of the meson wavefunctions that is not provided 
by the effective theory. 

In this respect soft-collinear effective theory applied to heavy-to-light decays at large 
recoil is rather different from HQET for B — ► form factors. In the latter case one ob- 
tains non-trivial form factors relations beyond leading order in A, because the heavy-quark 
flavour symmetry also relates the initial and final hadronic states. (The same symmetry 
is responsible for the absolute normalization of the Isgur-Wise form factors at zero recoil.) 
An analogous symmetry that relates light meson states is absent in soft-collinear effective 
theory. 



6 Conclusions 

The effective field theory approach to soft-collinear factorization initiated by |2|, |[] provides 
a new framework for formulating the factorization theorems of QCD for processes with en- 
ergetic, nearly massless particles. Beyond this it opens up the possibility to investigate 
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power corrections to such processes more systematically than in the diagrammatic frame- 
work. In this paper we have taken a first step in this direction by deriving the effective 
Lagrangian to second order in the expansion parameter A, and the corresponding power 
corrections to weak interaction quark currents that convert a heavy quark into energetic 
light partons. 

We have formulated the effective theory in a position space representation rather than 
the less familiar hybrid position-momentum space representation adopted in ||. Although 
the physical content of the two formulations should be identical, the choice of representation 
implies significant technical differences. In the hybrid representation the rapid variations of 
collinear fields are extracted and become field momentum labels, such that a collinear field 
£ p (x) has only a residual ultrasoft spatial variation. In position space collinear fields 
are characterized by their most rapid spatial variation in every direction with no possibility 
to extract superimposed slow variations. In this formulation the multipole expansion of 
the argument of ultrasoft fields accomplishes the equivalent of extracting the ultrasoft 
variation of collinear fields in the hybrid representation. A consequence of this is that the 
gauge symmetry of QCD is implemented differently in the two formulations of soft-collinear 
effective theory. 

The main results of this paper can be summarized as follows: 

- We have derived the soft-collinear effective Lagrangian to order A 2 including ultrasoft 
quark fields. We have also found that the tree-level Lagrangian is exact, since Lorentz 
invariance protects every operator (at any order in A) from acquiring a non-trivial 
coefficient function. 

- We have derived the effective heavy-to-light current relevant for exclusive B meson 
decays to order A 2 and investigated its properties under reparameterizations of the 
arbitrary vectors that define the momentum decomposition in the effective theory. 



An expression for the current at order A has been previously given in ||20|| . We differ 
from this result by an additional term that comes from the coupling of transverse 
collinear gluons to the heavy quark. 

We then investigated whether relations between the transition form factors for B meson de- 
cay into a light pseudoscalar or vector meson survive when one includes power corrections. 
Including terms of order A we find one relation among the three form factors for decays 
into pseudoscalar mesons and two relations among the seven form factors for decays into 
vector mesons. Including terms of order A 2 no relations survive. This is to be compared 
to the existence of two (five) form factor relations for pseudoscalar (vector) mesons in the 
infinite quark mass limit. 



Note added in proof. In the revised version of (|hep-ph/0205289| v2) Bauer, Pirjol 



and Stewart have also given the tree- level matching coefficient of the operator J 1; which 
appears in the expansion of the heavy-light current at order A, see (|120|) . Their result 
agrees with ours. 
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A Derivation of equations (|5t|) and fl6Q) 

In this appendix we provide the details of the manipulations that lead to the final form of 
the second-order Lagrangian, fl5"9|) and (|K]), for abelian gauge fields. We start from ([52]) 
and the explicit expressions for the various terms given below (|52|). 

To obtain (|58|) we use the equation of motion first for £, then for £. We have to keep 
the order A terms in the equation of motion and account for the fact that the Lagrangian 
changes after the first use of the equation of motion. After performing the steps that lead 
to ( |58"D one has to add to (|53|), (|54]) the terms 



AC 



(2) 



£ (gx±A us ) [in^D + iIp Lc 
1 



■ gA±us + gAi 



l 



AC 



(2) 



tj j _1_ LIS ' tj / _1_ LIS • j-^ 

m + D c m + D c 
Z(igx±A us )iip ±c W c q+ "h.c." 



ilp ±c , igx ± A v 



4- 



(147) 



from the equation of motion. Formally, using the equations of motion in ( |57D is equivalent 
to the change of variables £ — > (1 + igx±A ns ) £ in >Q and C^. 

(2) 

The terms in AC^ can be simplified by using the identities (valid only for abelian 
gauge fields) 



(x±A w ) (in-D) (x ± A n8 ) = - {(m_£>), (x ± A ns ) 2 } 



{x±A m ) iip. 



±c 



in + D, 



ijp. 



±c 



2 in+D, 



4us> {ixxA v 



-ilpxc (xj_A ns ) 

4-Lus 



in + D c 



in + D c 



in + D, 



in + D c 



(148) 



Using again the leading order equations of motion for £ and £, one obtains 



AC 



(2) 



_l_us 



in + D c 



(149) 
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which exactly cancels the corresponding term in (|53|). The remaining terms in >Q can be 
simplified by commutator relations and the equation of motion. The identity 



{in+D c y 



in + D c 2 



n_2 



(150) 



allows us to write 
1 



til?. 



_Lc 



in + D c 



{gn + A us )-^-iip^i= \Un-x) [{n_d) (gn + A us )} *±£. (151) 



This term combines with the first term in the third line of (^) into n l ^n u _F^. Then we 
use 



2 [( X± d) 



lJp±_ c + ilp ±c 



in + D c 



in + D c 



[(x ± d) 



(152) 



in + D, 



■tip. 



±c 



in + D c 



flU 



to rewrite the second line of (|53|) : 
1 



in + D c 



[{x±d) gA ±us ] + [(x±d) , 



s f ^ * M + * M ™k t « 



(153) 



With this all remaining terms involving ultrasoft gauge fields can again be expressed in 
terms of field strength tensors and (|59"D follows. The derivation of (|60"D , where one has to 
take account of ( |147|) , is straightforward. The application of the £ equation of motion at 
several steps in this calculation is equivalent to the field redefinition (|6lD . 



B Heavy- to- light form factors 

The form factors relevant to heavy-to-light B decays are defined as in 



(P(p')\^ 1 %\B(p)} = f + (q 2 



2 2 

p^+p 1 * 2— — tLq^ 



fo(q 



2 _ ™2 

2 T> 



2 m B m P 



(P(p')\^a^q u b\B(p)) 



ifr{q 2 ) 
mp + mp 



g 2 (p M + p ,f *) - (m B - mp) q 



2 \ 



(154) 
(155) 
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where is the B meson mass, m B the mass of the pseudoscalar meson and q = p — p'. 
The relevant form factors for B decays into vector mesons are defined as 



(V(p',e*)\^b\B(p)) = 2lV{ f ] e^el P ' p p a , 

rriB + mv y 



(156) 



£ ■ q 



(V(p',e*)\^Yl 5 b\B(p)) = 2m v A (q 2 )^q» + (m B + m v )A l (q 2 ) 



m B + my 



q* 



e ■ q 



(157) 



(V(p',e*)\i>a^q v b\B(p)) = 2r 1 (g 2 )e^£>X 



(158) 



(V(p',e*)\^ l5 q u b\B(p)) = H)T 2 (q 2 ) (m 2 B - m 2 v ) e*» - (e* ■ q) {jf + p'») 



+ (-i)T 3 (q 2 )(e*.q) 



TTln — m 



ip^+p'") 



B " V 



(159) 



where my (e) is the mass (polarization vector) of the vector meson and we use the sign 
convention e 0123 = —1. 

Calculating the various traces in ( |136| ) we can express the form factors ( 154j) to ( |1 59|) 
in terms of the sets fp and f± , f» introduced in (|137|) . Neglecting terms suppressed 
by a relative factor m 2 PV /E 2 ~ A 4 we obtain 



f + (q 2 ) = ff + \ E )-f P + -\E)-\(f p - + \E)-f p --\E) 
H<?) = ^^{f P ++ \E)-f P + -\E) + t r {f p - + \E)-f p --\E) 



m 



B 



13 



fAq 2 ) = T ^^{f P ++ \E) + f P + -\E)-f P ^\E)-f P ~-\E)} 
for decays into light pseudoscalars, and 

V(q 2 ) = r ^^{ft + \E) + f i r\E)-ft + \E)-f { r\E)}, 



m B 

Mf) = ?^{4 ++) (^)-4 + - ) (^) + 4M" +) (^)-4" ) (^) 

2m B my [ 11 11 m B x " " 



A 2 (q 2 ) 



m B — q 



m B (m B + m v ) 



[ff + \E) + f£~\E) + f[- +) (E) + f[-\E)} 



m B + m v j (++) 



- f^w - \ (fir + \E) - fi~-\E)) 

m B 



m B 



(160) 
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-t- 

m B 


Uq 2 ) 


= f[ ++ \E)~ 


T2(q 2 ) 


m 2 B -q 2 f 




T 3 (q 2 ) 






+f[ ++ \E) 



m 



B 



A ++ \E) - f[ + -\E) + 4 {ft + \E) - f[-\E)) 



T?\ , f(- + )( Z?\ , X( )/ 



ni 



B 



(161) 



for decays into light vector mesons. 
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